Predicting sport fish mercury contamination in heavily managed reservoirs: Implications for human and ecological health

Mercury (Hg) is a concerning contaminant due to its widespread distribution and tendency to accumulate to harmful concentrations in biota. We used a machine learning approach called random forest (RF) to test for different predictors of Hg concentrations in three species of Colorado reservoir sport fish. The RF approach indicated that the best predictors of 864 mm northern pike (Esox lucius) Hg concentrations were covariates related to salmonid stocking in each study system, while system-specific metrics related to productivity and forage base were the best predictors of Hg concentrations of 381 mm smallmouth bass (Micropterus dolomieu), and walleye (Sander vitreus). Protecting human and ecological health from Hg contamination requires an understanding of fish Hg concentrations and variability across the landscape and through time. The RF approach could be applied to identify potential areas/systems of concern, and predict whether sport fish Hg concentrations may change as a result of a variety of factors to help prioritize, focus, and streamline monitoring efforts to effectively and efficiently inform human and ecological health.


Introduction
Mercury (Hg) is a concerning contaminant due to its widespread distribution and tendency to bioaccumulate in organisms to concentrations that negatively influence the health of humans, other organisms, and ecosystems worldwide [1][2][3]. Mercury enters the landscape through a variety of natural and anthropogenic pathways including volcanoes, forest fires, erosion, fossil fuel burning, waste incineration, mining operations, and cement production [4]. Despite health concerns related to Hg contamination in fish, documented benefits of fish consumption may outweigh the risks [2,5,6], and healthy fish consumption has been encouraged [7]. Thus, understanding fish Hg concentrations and monitoring programs are necessary to protect human health related to fish consumption.
Sport fish Hg concentrations are generally a preferred indicator of Hg contamination in freshwater systems versus other aquatic organisms. Consumption of contaminated sport fish (with Hg incorporated in their muscle tissue) is the major method of transfer of aquatic Hg to humans [1,2]. However, fish Hg concentrations can be inherently variable due to several factors, and changes in contributions of Hg to the environment do not necessarily translate directly to changes in fish Hg concentrations because of complex biogeochemical interactions [8,9]. Thus, Hg deposition information alone cannot be used to infer fish Hg concentrations across the landscape, and extensive, large-scale monitoring programs can be cost prohibitive.
Resource limitations within agencies, variability in fish Hg concentrations, and complex environmental interactions, make predictive model development (to identify priority areas or species of concern for testing) challenging. Approaches incorporating information about mechanisms influencing fish Hg concentrations and multiple data sources (including correlated predictors) could help overcome some of the obstacles to characterizing dynamic uncertainty in fish Hg concentrations. Investigators have studied a variety of factors associated with Hg contaminated sport fish in North America [10][11][12]. Several patterns have emerged on the larger landscape. For example, when deposited, Hg can be taken up quickly by water and the organisms within it [13]. Because of this direct water-atmosphere interface, larger systems (higher surface area) receive more overall Hg deposition relative to smaller systems. In addition, systems in close proximity (<100 km) to "large" sources of Hg (similar to those found in the eastern and central United States) can experience increased Hg deposition [14]. Wentz et al. [15] found that Hg inputs to stream fishes were highest in urban versus rural areas. It follows that systems located near Hg sources (i.e., population centers or urban areas versus rural areas) could experience higher levels of Hg inputs and exposure.
Fish species, size, and trophic position are important factors influencing sport fish Hg concentrations [16][17][18]. Specifically, large, piscivorous sport fish with elevated trophic positions tend to have higher Hg concentrations than small, omnivorous or planktivorous fish feeding at lower trophic levels. Therefore, changes in food web structure (especially the forage base of sport fish) can alter sport fish Hg concentrations by changing fish trophic position, growth, and diet [19][20][21]. Aquatic food web structure and species interactions are often influenced by fisheries management practices (e.g., harvest regulations, fish stocking). Thus, these potential sources of variability are important to consider.
In the arid portion of western North America, artificial reservoirs are the predominant lacustrine systems on the landscape, and tend to be highly managed with respect to their biotic and abiotic characteristics. These characteristics can influence Hg dynamics [22], and several studies have found that reservoirs have fish and other biota that can be particularly high in Hg concentrations in general [23][24][25]. Reservoirs can experience changes in water levels, and water level fluctuation has been associated with elevated Hg concentrations in sport fish [26,27]. These studies suggest that fish Hg levels are influenced by water level fluctuation through the rewetting and perturbation of dry soils which are relatively rich in sulfate, stimulating population growth of sulfate-reducing bacteria, Hg methylation rates, and Hg uptake by biota. Analogous patterns of elevated Hg concentrations in biota (zooplankton and fish) in systems experiencing water level fluctuation have also been found [28].
Nutrient inputs also influence system productivity and subsequent Hg concentrations in sport fish. Increased productivity can reduce Hg concentrations in organisms through "bloom "dilution [29,30] or "growth dilution" [31]. When high-quality diet items are available to organisms, they typically display high rates of somatic growth paired with lower consumption rates. This reduces Hg concentration in prey and overall intake by predators. In the case of fish, it has been shown through experimental nutrient additions and observational studies that fish from systems with higher nutrient inputs have relatively low Hg concentrations [32][33][34]. Thus, increased nutrient inputs to reservoirs can reduce Hg bioaccumulation in sport fish within them. At the same time, Hg is methylated (becoming bioavailable) under anoxic conditions by bacteria as a byproduct of their energy sequestration pathway [35,36]. High nutrient inputs can stimulate primary production and decay creating anoxic conditions conducive to Hg methylation [37,38]. Hypoxic conditions have been associated with high concentrations of Hg in water, zooplankton, and fish [11,39].
Several additional factors can act as drivers that influence fish Hg concentrations depending on environmental or other conditions. For example, water temperature can influence growing season length (and subsequently growth dilution), community composition, and Hg methylation rates. All else being equal, fish that grow faster can be "diluted" in Hg relative to others growing more slowly [40]. Thus, more growth and longer growing seasons can result in lower fish Hg concentrations. However, warmer water may also stimulate microbial activity and Hg methylation, making Hg more available for uptake by biota. Therefore, elevation and latitude could influence sport fish Hg concentrations positively or negatively as they relate to differences in system-specific temperature regimes.
We describe a machine learning approach to identify predictors of Hg concentrations in three different sport fish species in Colorado, USA reservoirs. We use data compiled from sources throughout Colorado to test for the relative importance of various predictors of sport fish Hg concentrations related to the factors described above. We used available system-specific predictors as indicators of urbanization, food web structure, productivity, and Hg methylation processes. Further, we used our observations combined with those of others to compare rates of sport fish Hg change that might be expected in response to a variety of factors primarily focusing on empirical changes in food webs and hypothetical changes in Hg emissions. Our comparisons provide a basis to inform decisions about the magnitude and relevant spatial and temporal scales for developing appropriate fish consumption advice.

Data collection and compilation
Sport fish Hg concentrations (2004 to 2013) available from Colorado Department of Public Health and Environment (CDPHE) were used as response variables for three fish speciesnorthern pike (Esox lucius), smallmouth bass (Micropterus dolomieu), and walleye (Sander vitreus)-across 32 reservoirs (Fig 1). For each system-species combination, simple linear regression was used to estimate fish Hg concentrations (and upper 95% confidence intervals when regressions were based on five or more values) at lengths that were related to the most commonly applied length-based harvest regulations (i.e., minimum length limits) in Colorado (864 mm for northern pike and 381 mm for smallmouth bass and walleye). We expected sport fish Hg concentrations to increase, decrease, or remain the same with increasing lengths depending on the characteristics of a given reservoir. Thus, estimates of fish Hg concentrations at their respective, relevant lengths were used whether the regression slope was positive, negative, or not significantly different than zero. To develop the regression analyses and determine Hg at the specified lengths, data from individual fish and composite samples were used. Composite samples were represented by a single datum as the mean length of the group of fish tested and the resulting Hg concentration. Individual fish lengths in composite samples were within 5% of the measured mean for all fish included in any given composite sample. Measurements of Hg concentration that fell below detection limits (variable across years) were censored by assigning a value equal to half the detection limit [41]. Regression analyses for single system-species combinations were never based solely on data that fell below a Hg concentration detection limit of 0.3 ppm because it was felt the assigned value (0.15 ppm) could mask a relatively large amount of variability (potential empirical range from 0 to 0.3). In 26 of 730 cases (less than one half percent of the data available), censored values (0.15 ppm) for the detection limit of 0.3 ppm were applied. In these cases there were additional data available to corroborate/support the censored data. Since censored values were similar to uncensored information, results were not sensitive to the inclusion or removal of these 26 values in regression analyses, or overall findings. When less than five individual or composite samples were available for a given system-species combination, those estimates were only retained if data from other species within the same system, or data from the same system but outside the 2004-2013 timeframe, corroborated the estimate. Samples were collected from different systems in different years and with different methodologies (i.e., sample sizes, composites versus individuals, species analyzed) and these are described in Table 1.
We selected a regression approach to maximize the amount of data that could be included in these analyses and to represent a length that is informative to the health of anglers and their families that harvest sport fish in Colorado. Our approach generally limited extrapolation outside Hg concentrations estimated for 864 mm northern pike and 381 mm smallmouth bass and walleye. However, extrapolation (using data from fish that did not encompass the length where point estimates were made) was necessary for some (described below) of the 54 reser-   mm) these combinations included Rifle Gap (461-525 mm), Stagecoach (388-457mm), and Totten (520-635 mm). These system-species combinations were retained because additional data (after 2013 or from other fish species) corroborated point estimates, and including these system-species combinations allowed us to produce predictions there. We note that this was done for prediction purposes to inform fish Hg testing, and would not be appropriate to develop fish consumption advisories to protect human health. We used data from 32 disparate reservoirs across Colorado with a variety of municipal, agricultural, and recreational uses. Predictors of the response variables (fish Hg concentrations at the specified lengths) were selected based on previous research, ease of data collection, and data availability. We primarily used data that were accessible to us, and also relatively easily collected or calculated metrics, so others could use analogous measurements across a variety of systems, potentially with data already available. We acknowledge that some of these metrics are oversimplified, but complex community, atmospheric and watershed-level modeling to capture Hg inputs and cycling were outside the scope of this study. Thus, metrics used for predictions should be selected with forethought, and potential confounding factors should be considered. We had a priori hypotheses associated with each predictor variable, and for each system, predictors included a group of 17 indices developed from a compilation of data from 2000 to 2014. These covariates included: 1) system surface area (Area; ha), 2) shortest linear distance to Colorado Interstate 25 (I-25 Data compiled for these analyses were collected intermittently by different agencies and personnel from 2000-2014 depending on a variety of factors including research or monitoring objectives, resource availability, and sampling conditions. We used these data to develop system-specific indices to characterize covariates and test for their relative importance for predicting empirical sport fish Hg concentrations. When possible, we limited data use from 2003- Each study reservoir is listed and sample sizes for northern pike (NPK), smallmouth bass (SMB), and walleye (WAL) are provided along with the year or years data were collected. Individuals that were tested for Hg concentration are indicated with a sample size followed by the letter corresponding to the first letter of the species represented. Composite samples are indicated with a sample size followed by the letter corresponding to the first letter of the species represented (N, S, and W for northern pike, smallmouth bass, and walleye, respectively) and a parenthetical range of individuals per composite sample. In systems where empirical data were not available for a given species, that system was excluded from that species analysis. Data that were censored are indicated with a sample size followed by the letter corresponding to first letter of the species represented and a parenthetical detection limit relevant to those data. Estimates of fish Hg concentrations (ppb) from empirical data (Obs) are provided following the letter corresponding to the first letter of the species represented as well as the upper 95% confidence interval (parenthetically, following the Hg concentration) when five or more samples were available to calculate error for a given system-species combination. Where error could not be calculated, only the estimated Hg concentration (ppb) is provided parenthetically without the upper 95% confidence interval. https://doi.org/10.1371/journal.pone.0285890.t001 2013 to develop the most relevant covariates for the 2004-2013 Hg concentration data for Colorado sport fish. We did this to ensure that sport fish Hg concentrations were temporally reflective of the system-specific conditions characterized by the indices/covariates. However, contemporaneous data were not always available, and supplemental data from 2000-2002 and 2014 were necessary to characterize some systems (see Table 2 for covariate descriptions). Thus, our indices/covariates were developed to inform the relative importance of past drivers of sport fish Hg concentrations, and design future monitoring efforts, but not to develop fish consumption advice directly. Data from multiple system-specific surveys and sampling events were available from CPW   were censored by assigning a value equal to half the detection limit [41]. Values/categorizations of these metrics are provided for clarity and comparison (Table 3). Data were acquired from public sources, and fish Hg concentration data are available for download from CDPHE (https://cdphe.colorado.gov/water-quality-records-center-and-requests). The predictors of fish Hg concentrations tested for each reservoir are provided in Table 3.

Data analyses
We applied the machine learning technique called random forest (RF) to statistically test for indices/covariates as predictors of sport fish Hg concentrations within the dataset and to test for the relative importance and effectiveness of the predictors for characterizing sport fish Hg concentrations [42]. The RF approach is ideal because it is based on a form of cross-validation, is entirely nonlinear, accounts for obscure interactions, and simultaneously evaluates multicollinear variables [42,43]. Although many other approaches prevent the inclusion of multicollinear variables, the RF approach can incorporate all readily available data for prediction. For example, TP and Chl a are often highly correlated, and with most approaches, one variable or the other would need to be excluded from any single predictive model. Using the RF approach allowed us to incorporate all variables of interest despite potential multicollinearity. The RF approach falls under a much larger class of machine learning methods where prediction is of primary interest. These methods rely on "bagging," which refers to the method of fitting regression trees to numerous bootstrapped versions of the training data (i.e., observations sampled with replacement from the larger data set). The resulting trees are then averaged to yield a predictor with low bias and variance. To perform our RF analyses, we used the 'ran-domForest' package 4.6-6 [44] in R [45]. In the RF implementation described here, 2,000 regression trees were used to calculate accuracies and error rates for each observation using out-of-bag predictions (i.e., predicting data that were withheld from each tree). The predictions of data that were not used to evaluate fit can be considered as a form of cross-validation. Variable importance can then be assessed by comparing the increase in; 1) mean squared prediction error, and 2) node purity associated with each individual covariate. We refer the interested reader to several works for more information on these variable importance metrics [44][45][46].
For each fish species, a full model (all predictors) was used to test for predictive power and variable importance. Although we selected metrics that were relatively common and available, we were aware that some of the predictors used may not be available. To determine how well this approach performed when only the best predictors were used, we removed most predictors from the full models and repeated the RF statistical analyses on a condensed model for each species. Predictors were retained if they were found to be important as measured by the increase in mean square prediction error and node purity. These selections were made somewhat arbitrarily for demonstration purposes, however, a priori hypotheses, and support from other research were considered along with the predictor importance metrics when full models were condensed to the most important three or four variables. For each system-species combination, we compared the estimated Hg concentrations from empirical data to the Hg concentrations estimated from the full RF model. The difference of these (empirical estimate minus RF prediction) was regressed as a function of the empirically estimated Hg concentrations to test whether this differed from a one-to-one relationship to evaluate the presence of any bias. In general, we expected that the best predictors of Hg concentrations in 864 mm northern pike would be those related to stocked salmonids as forage while the best predictors of Hg concentrations in 381 mm smallmouth bass and walleye would be more dependent on system productivity metrics.

Prediction application
The food web structure in Elkhead Reservoir (Colorado, USA) changed significantly when the stocking of salmonids (up to 3,400 kg of~250 mm rainbow trout stocked in a single year; 2010, with a mean of about 1,500 kg of salmonids stocked each year) was discontinued in the fall of 2011. Based on previous whole-system experimental work [20] we expected northern pike Hg concentrations in this reservoir to increase from this management action. For comparison, northern pike Hg concentrations were available from 2005 and 2007, prior to this management shift when significant masses of catchable (~225 mm) rainbow trout were being stocked. We used the full and condensed models developed for northern pike to predict what 864 mm northern pike Hg concentrations would be under these new conditions (setting all salmonid stocking to zero). The output provided a value for comparison to empirical northern pike Hg concentration data that were collected late in 2013 (following the cessation of stocking), that was not included as training data for the RF approach. This provided information from the same system with and without salmonid stocking (a factor expected to have predictive importance for large northern pike) to evaluate model performance, recognizing that the food web structure may or may not have been in equilibrium following this shift.

Random forest modeling
Salmonid stocking metrics best explained 864 mm northern pike Hg concentrations. Thirteen systems were included in the 864 mm northern pike RF analysis. The Hg concentration in northern pike was also predicted from Elkhead Reservoir experiencing no salmonid stocking. The full model explained about 25% of the variability in northern pike Hg concentrations with all three salmonid stocking metrics being the most important predictors of northern pike Hg concentrations (Fig 2). To reduce the number of predictive variables in the model, only the top three predictors related to stocking were included in a second model. This more parsimonious, condensed model (northern pike Hg concentration as a function of Salm mass, RBT kg/ha, and RBT) explained about 50% of the variability in the data. There was a consistent bias from a one-to-one line (slope; 0.33, F = 28.79, p-value < 0.01) where predicted values were higher than those observed at relatively low concentrations and lower than those observed at relatively high concentrations (Fig 3).
The best predictors of smallmouth bass Hg concentrations were metrics related to growing season, productivity, and forage base. Fifteen systems were included in the smallmouth bass RF analysis. The full model explained about 55% of the variability in smallmouth bass Hg concentrations (Fig 4). Although I-25 appeared as the fourth most important predictor in both variable importance indicators, the relationship was contrary to expectations; smallmouth bass Hg concentrations went up with increasing distance from Colorado Interstate 25. Thus, this predictor was removed from further analyses and the top three predictors (Elev, GSD, and Chl a) were used in the condensed model. The more parsimonious, condensed model (smallmouth bass Hg concentration as a function of Elev, GSD, and Chl a) explained about 70% of the variability in the data. There was a consistent bias from a one-to-one line (slope; 0.25, F = 17.78, pvalue < 0.01) where predicted values were higher than those observed at relatively low concentrations and lower than those observed at relatively high concentrations (Fig 5).
Walleye Hg concentrations were best predicted by metrics of system productivity and growing season. Twenty-six systems were included in the walleye RF analysis. The full model explained about 10% of the variability in walleye Hg concentrations (Fig 6). Thus, the top four predictors (Chl a, Elev, Secc, and TP) were used in the condensed model. The more parsimonious, condensed model (walleye Hg concentration as a function of Chl a, Elev, Secc, and TP) explained about 17.5% of the variability in walleye Hg concentrations. There was a consistent bias from a one-to-one line (slope; 0.35, F = 54.17, p-value < 0.01) where predicted values were higher than those observed at relatively low concentrations and lower than those observed at relatively high concentrations (Fig 7).

Prediction application
The linear regression-based empirical Hg concentration of 864 mm northern pike in Elkhead Reservoir based on 2005-2007 data (during salmonid stocking) was 0.02 ppm. When the full model was applied to predict Hg concentration of 864 mm northern pike in Elkhead Reservoir without stocking (setting Salm mass, RBT mass/area, and RBT to zero in the model), the result was 0.55 ppm. The result from the condensed model (including RBT, Salm mass, and RBT mass/area as predictors) when the same was done (stocking set to zero) was 0.93 ppm. The observed 2013 (salmonid stocking ceased) Hg concentration of 864 mm northern pike in Elkhead Reservoir was 0.57 ppm.

Random forest modeling
Our expectations based on previous work in Colorado reservoirs (20,21,47,48) regarding the factors important for predicting sport fish Hg concentrations were supported by the RF  Table 2. analysis. The most important predictors (based on the RF variable importance indicators) of 864 mm northern pike Hg concentrations in Colorado were related to salmonid stocking, presumably because of biomass dilution from the consumption of stocked fish that are relatively high in energy and low in Hg content [20,47,48]. We expected that the most important predictors of 381 mm smallmouth bass and walleye Hg concentrations would be related to biomass dilution as well, but driven by system productivity rather than fish stocking. This is likely because 381 mm sport fish are not generally capable of consuming large (> 200 mm) stocked salmonids [49] that represent the majority of the biomass of fish stocked in Colorado. Thus, indicators related to system productivity and growing season (i.e., Elev, Chl a, Secc, TP) were found to be relatively important predictors of smallmouth bass and walleye Hg concentrations compared to stocking indices.
Four variables that were included as predictors of sport fish Hg concentrations were found to be relatively unimportant in the RF analysis. The indicators of proximity to emission sources and deposition potential (Area, I-25, Pop dens) were relatively unimportant as was the indicator of water level fluctuation hypothesized to be positively correlated to Hg methylation [26,27]. Instead, we found that productivity (internal and external sources) in the form of nutrients and forage were better predictors of Hg concentrations in sport fish in this set of systems relative to indicators of Hg deposition and methylation within our dataset. Although studies have shown that Hg sources and deposition may be elevated in and around urban areas, that does not necessarily translate to elevated Hg concentrations in fish within those areas [15,50]. These observations (fish Hg concentrations decoupled from deposition/sources to some degree) could explain some of our findings related to urbanized areas in Colorado. In Colorado, system elevation is related to several other factors that could potentially increase and decrease fish Hg concentrations. As elevation increases agricultural land-uses decline along with nutrient inputs, precipitation (linked to Hg deposition) increases, growing season declines, and in-lake forage species (e.g., gizzard shad with relatively high energy and low Hg content) are less abundant, which could all lead to elevated sport fish Hg concentrations. However, with increasing elevation, food sources tend to have relatively low trophic positions, food chains tend to be relatively short, and there are generally fewer local sources of Hg emissions, which could all lead to lower sport fish Hg concentrations. In this case, factors  Table 2. https://doi.org/10.1371/journal.pone.0285890.g004

PLOS ONE
related to elevation that led to increased fish Hg concentrations appeared to have a stronger influence than those leading to decreased fish Hg concentrations. The precise mechanism or combination of mechanisms causing this pattern are unknown, but the other three most important predictors of fish Hg concentrations were indicators of system productivity.
An RF approach analogous to what is described here is appropriate to inform and streamline monitoring programs, but not to obtain estimated Hg concentrations as an alternative to direct testing. This technique could be used to identify areas where sport fish Hg concentration shifts may be occurring. Appropriate systems and species could then be prioritized for Hg testing so consumption advisories could be updated accordingly to maximize the benefits of fish consumption and better protect anglers and their families from potential health risks. For example, in response to the management actions and subsequent food web shifts observed in Elkhead Reservoir, testing was prioritized and CDPHE has adjusted fish consumption advice and monitoring efforts there to be more protective of human health.

Prediction application
The performance of the full and condensed RF models was encouraging. Although there was a range of variability described by the full and condensed models, model results were reasonable despite low sample size, and bias appeared systematic (predicted values were consistently higher than those observed at relatively low concentrations and lower than those observed at relatively high concentrations). The application of the full model to Elkhead Reservoir to predict northern pike Hg concentrations following 2011 (when salmonid stocking was discontinued), resulted in a value within 5% of empirical observations from 2013, while the application of the condensed (three stocking variables) model overestimated northern pike Hg concentrations. Thus, the full model set performed well for making predictions in Elkhead Reservoir specifically, but the condensed model set was more parsimonious and explained more variation in the data across all of the reservoirs considered during the northern pike analysis. It is likely that within Elkhead Reservoir, other predictors included in the full model (in addition to salmonid stocking), may be important for predicting northern pike Hg concentrations, though  Table 2.
https://doi.org/10.1371/journal.pone.0285890.g006 further investigation at the system-specific level would be necessary to identify mechanism(s). This finding highlights the potential value of the RF approach and including a variety of factors during predictive analyses.

Conclusions
Many reservoirs in western North America are heavily managed, highly fluctuating, and serve multiple uses which can influence their Hg dynamics [22]. Biota in reservoirs in particular have been found to have elevated Hg concentrations [23][24][25]. These types of heavily managed and dynamic systems are analogous to those considered here, and Colorado offered a wide range of predictors like high and low elevation and little or no salmonid stocking to heavy salmonid stocking. Thus, the patterns observed in this study may be present in other regions where comparable lake and reservoir management approaches (e.g., salmonid stocking, irrigation, recreation) are being applied. The RF approach described here could identify important predictors of fish Hg concentrations at other locations and scales across the landscape to inform monitoring and advisory programs in place to protect human health.
Based on our observations, sport fish Hg concentrations varied widely across the landscape, driven by multiple, species-specific factors. Thus, simplistic yet flexible approaches to rapidly identify areas or species of concern, and communicating findings to anglers and their families (e.g., the case of the Elkhead Reservoir northern pike population), are likely the most effective means for addressing human health concerns from Hg contamination in the short-term. Based on available atmospheric information [14,51,52], Hg deposition may increase in some areas (including areas of Colorado) in the future. These potential changes coupled with the observed variation in sport fish Hg concentrations, make the RF approach a useful, predictive tool to prioritize and streamline Hg monitoring efforts. Control of anthropomorphic Hg emissions represents an important component of the ubiquitous issue of Hg contamination in the environment. This is recognized by groups like the Minamata Convention on Mercury, and likely represents one of the few long-term solutions to the problem of Hg significantly exceeding natural concentrations in ecosystems and biota. In the short-term, predictive tools like RF that incorporate food web structures, and their dynamics, are important to understand and provide context for Hg bioaccumulation in sport fish and changes in Hg emission, deposition, and cycling.